A bottleneck model for bidirectional transport controlled by fluc- 
tuations 
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Abstract - We introduce a new model to study the oscillations of opposite flows sharing a 
common bottleneck and moving on two Totally Asymmetric Simple Exclusion Process (TASEP) 
lanes. We provide a theoretical analysis of the phase diagram, valid when the flow in the bottleneck 
is dominated by local stationary states. In particular, we predict and find an inhomogeneous high 
density phase, with a striped spatio-temporal structure. At the same time, our results also show 
that some other features of the model cannot be explained by the stationarity hypothesis and 
require consideration of the transients in the bottleneck at each reversal of the flow. In particular, 
we show that for short bottlenecks, the capacity of the system is at least as high as for uni- 
directional flow, in spite of having to empty the bottleneck at each reversal - a feature that can be 
explained only by efficient transients. Looking at more sensitive quantities like the distribution of 
flipping times, we show that, in most regimes, the bottleneck is driven by rare fluctuations and 
descriptions beyond the stationary state are required. 
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The understanding of the macroscopic behavior of com- 
plex systems out of equilibrium is one of the main chal- 
lenges in modern statistical mechanics. A common feature 
of many non-equilibrium systems is the presence of a cur- 
rent in their stationary state, in contrast to equilibrium. 
A general framework describing these systems is still lack- 
ing, though the importance of current large deviations and 
their link with fluctuations theorems has been emphasized 

The absence of such general framework motivated the 
study of many oversimplified microscopical models, among 
which exclusion processes have been reference systems be- 
cause they allow for extremely precise numerical results 
and exact analytical solutions in some cases. Exclusion 
processes are simple models defined on a discrete -usually 
one-dimensional- lattice, on which particles hop from site 
to site. Their role as a reference system is of particular 
importance for exact calculations of large deviations, the 
non-equilibrium counterpart of the free energy (2j|3]. At 
the same time, exclusion processes are a flexible tool to 
model various physical systems. Indeed, they capture the 



correct collective behavior of systems with very different 
length scales from social individuals, such as pedestrians 
4 , vehicles 5 , ants [4[[6] , to molecular systems as molec- 
ular motors 3 [i] , and microscopical ones as quantum dots 
ff. Furthermore exclusion processes are closely related 8] 
to growth phenomena described by the celebrated Kardar- 
Parisi- Zhang equation [9], and the universality of these 
phenomena has been confirmed in recent experiments in 
liquid crystals 10 . 

For the asymmetric simple exclusion process (ASEP) 
the full current probability distribution has been calcu- 
lated for different boundary conditions 11,12 . However, 



while the ASEP captures the proper behavior of particles 
moving on a single lane, it is obviously not appropriate to 
describe situations in which the particles (such as individ- 
uals or cars) move in few intersecting lanes. To describe 
such systems several one-dimensional exclusion processes 
must be coupled and the task of determining the current 
distributions becomes harder. Models for interacting par- 
allel lanes (to describe e.g. the traffic on highways) have 
been introduced and studied [l3 . In particular, in the 
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Fig. 1: Schematic representation of the model. 



so-called bridge models, two lanes share a finite number 
of sites (the 'bridge') [14^. Such systems have attracted a 
large interest due to the symmetry breaking that occurs 
in most cases MjjlS]. 

In this letter, we introduce a new model which couples 
two totally asymmetric simple exclusion process (TASEP) 
lanes with oppositely directed flows sharing a common 
bottleneck (see Fig. [T]). In contrast to the bridge models, 
exchange of oppositely moving particles is not possible in- 
side the bottleneck, because we impose that only particles 
going in one direction can go through the bottleneck at 
a given time. Therefore, particles going in the opposite 
direction have to wait until the whole bottleneck is empty 
before being allowed to go through. Our model can be seen 
as a representation of opposite pedestrian flows crossing 
e.g. at a door. The model is also relevant for other systems 
such as multiphase flows, or bidirectional molecular traffic 
across nuclear pores 
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We study this model with a combination of analytical 
and numerical techniques. The property that makes this 
model qualitatively different from previous ones is that the 
condition for reversing the flow inside the bottleneck (i.e. 
for having an empty bottleneck) is a rare event as soon as 
the bottleneck exceeds a few sites. Indeed, the dynamics 
of this system is driven by rare fluctuations, and is non 
trivial. For example, there exists a regime in which stop 
and go waves invade the whole system. Also, a counter- 
intuitive feature of the model is that it can sustain in the 
bottleneck a current higher than the maximal current of a 
single lane system of infinite size. 

The model. — We define now the model more pre- 
cisely. Particles move on two parallel tracks, modeled as 
two TASEP lanes. Both lanes share the same bottleneck 
of length Ni). We call (resp. '-') the particles moving 
from left to right (right to left), represented as red (blue) 
particles in Fig. [l] We consider explicitly only the lanes 
of incoming particles (of length L). Thus outgoing lanes 
are shadowed in Fig. [l] Random sequential update is con- 
sidered. Particles enter the lane with rate a, hop forward 
with rate p if the next site is empty, and leave the bot- 
tleneck with rate p. A particle enters the bottleneck with 
rate po only if (i) there is no particle of the other species 
inside the bottleneck and (ii) the first site of the bottle- 
neck is empty. We assume p = po = 1, leaving only a and 
/3 as the model parameters, in addition to the bottleneck 
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Fig. 2: Phase diagram of the bottleneck model, predicted by 
the phenomenological approach. 



and lanes lengths and L. 

Phase diagram. — By using the stationary proper- 
ties of the TASEP, we infer the phase diagram (given in 
Fig. |2| of the system with bottleneck from simple phe- 
nomenological arguments. In the following, currents refer 
in general to the current of one given species of particles 
(and not the total current). 

For low a values, we expect a free flow (FF) phase driven 
by the entrance rate. For each species, the corresponding 
particle current and density are 
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For large a and small /3, we expect an exit driven jammed 
(JJ) phase. When a species goes through the bottleneck 
for a long enough time, the system should become equiv- 
alent to a one-lane TASEP in stationary state, and the 
current should be /3(1 — /3). As the use of the bottleneck is 
shared between the two types of particles, each species can 
go through only half of the time. Neglecting the effects of 
the transients, on average the current is 



(2) 



We shall see that the particle distribution in incoming 
lanes is not homogeneous, and that, in contrast to cur- 
rent, density cannot be obtained from this mean field like 
approach. 

The FF/JJ boundary is given by the solution of Jin = 
Jouti satisfying ac{P) when /3 ^ 0, i.e. 



a, 



:(/^) = ^-^Vl- 2/3(1-/3). 



(3) 



The FF phase corresponds to a < ac{l3). 
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Fig. 3: Average current as a function of a, for different values 
of f3. Symbols indicate numerical results obtained from Monte 
Carlo simulations for L = 50 and Nh — 4. The solid line 
gives the theoretical predictions for /S > 1/2, i.e. both for the 
free flow and the maximal current phase. Dashed lines give 
predictions for /3 < 1/2 for the jammed phase with a > ac{/3). 

When a > 1/2 and /3 > 1/2, the single lane TASEP is 
in the maximal current (MC) phase, with a current equal 
to 1/4. Again, if we neglect transients in the bottleneck, 
the average current must be half of this value: 

Jmc = 1/8. (4) 

The MC/FF boundary, obtained from Jin = JmCi is given 
by ac = (I - \/l72)/2 ^ 0.146, while for the MC/JJ 
boundary, J out = Jmc provides /3c = ^. 

Monte Carlo simulations confirm these mean field pre- 
dictions both for current and density in the FF phase [Eq. 
([T])], and for the current in the J J phase [Eq. ([2|] for 
values of P not close to the MC boundary (see Fig. [3|. 
By contrast, the current in the MC phase is systemati- 
cally underestimated by this phenomenological approach. 
However, as seen in Fig. [4j the maximal current value Q 
is recovered for long bottlenecks, for which a stationary 
state can be established in the bottleneck. For shorter 
bottlenecks, currents higher than the stationary maximal 
current 1/4 can be obtained. 

Interestingly, with our bi-directional model, we can even 
reach the capacity of a finite-size system of size 11 
which would be fed by two incoming lanes, each trans- 
porting the maximal current and injecting particles on the 
same end of the bottleneck (in this case, flow would be uni- 
directional in the bottleneck and no flux reversal would be 
needed). 

Hence the counterintuitive result according to which it is 
at least as efficient to have particles coming from both ends 
of the bottleneck than having all these particles coming 
from the same end. This high capacity can be obtained in 
spite of the fact that bidirectional traffic requires to empty 
the bottleneck at each flux reversal, a limitation which is 
compensated by very efficient transient states. Indeed, at 
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Fig. 4: Average current against /3 for different Nb, and L = 50, 
a = 0.6. The dashed line gives the theoretical predictions 
for a = 0.6, both for the jammed and MC phase. For large 
bottlenecks the current saturates to Jmc = 1/8. 

each reversal of the current in the bottleneck, particles 
enter an empty system. Their motion is not hindered by 
predecessors, and thus high fluxes can be achieved. 

Spatio-temporal structures. — To characterize the 
dynamics of the system, we use a domain wall approach 
[T7[|l8) describing the system at a mesoscopic scale. We 
consider only the FF and J J phases, since the domain 
wall approach is not appropriate to describe the MC phase 
(which has long-range correlations). In this approach, we 
neglect transients inside the bottleneck, and assume that 
currents and densities are given by stationary expressions. 

First we describe the FF phase, and take the point of 
view of particles. When the bottleneck is closed, 
particles accumulate in front of the bottleneck, forming a 
queue of density 1. The upstream end of the queue can 
be seen as a discontinuity (or a wall) separating the queue 
from the bulk. A new wall is created when the bottleneck 
opens. The queuing particles feed the bottleneck with a 
high effective injection rate, and a high density domain 
is installed in the bottleneck. Ignoring transients, this 
jammed domain imposed by the exit has density 1 — P 
and current /3(1 — /3). Due to mass conservation, the wall 
between the jammed domain and the queue moves back- 
wards with velocity 1— /3 until the whole queue is dissolved. 
Then, at the separation between the bulk FF phase and 
the jammed domain localized at the exit, a new wall forms 
and moves forward with velocity V = ^ 
order to understand which phase invades the bulk, we have 
to determine whether the FF domain will reach the en- 
trance of the bottleneck before it closes again. This is the 
case if 

a(l - a) < - /?), (5) 

and assuming that the bottleneck is open and closed for 
the same period of time r (which is true on average). Then 
the queue and the jammed domain stay localized near the 
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Fig. 5: Spatio-temporal plot in the jammed phase {a = 0.6, 
/3 = 0.2, L = 100, Nb = 2). Both incoming TASEP lanes are 
shown, with the bottleneck in the middle. Particles (empty 
sites) are represented as black (white) spots. 

bottleneck, and the bulk FF phase can be sustained. If 
condition ([5| is not fulfilled, the bulk FF phase is slowly 
invaded by the queue, and cannot survive over long times. 
Note that Eq. ([5| is indeed identical to condition (|3| for 
the FF/JJ boundary. 

Now we apply the same coarse-grained approach to the 
JJ phase to show that a homogeneous bulk density is not 
possible. Indeed, pbuik should satisfy pbuik{'^ - Pbuik) = 
J out' However, when the bottleneck is closed, a queue 
of density 1 is formed. The wall between this queue and 
the bulk density moves backwards with velocity phuik- In- 
stead, when the bottleneck is open, the wall separating the 
exit driven jammed domain and the queue moves back- 
wards only with velocity 1 — j3. Thus, it cannot catch up 
with the previous wall, and the queue can never be en- 
tirely dissolved. When the bottleneck closes again, a new 
queue of density 1 is formed, whose rear end also moves 
with velocity 1 — (3 inside the exit driven phase. Then, 
regions of queuing particles with density 1 alternate with 
regions of density 1 — P and current /3(1 — /3), resulting 
in an overall striped jammed phase. This is confirmed by 
spatio-temporal plot in Fig. [5] obtained from simulations. 
Averaging over the stripes gives an actual bulk density 

Pont = l-^, (6) 

also confirmed by the Monte Carlo simulations. 

Distribution of the oscillation periods. — Until 
now, we considered the average value of the oscillation 
period r, which is a fluctuating variable. To perform a 
quantitative analysis, we define T as the time during which 
the bottleneck is occupied by at least one particle of the 
type under consideration. After T, bottleneck is empty 
and a particle of either the same or opposite type can 
enter. Thus T is not identical to r but strongly related 



Fig. 6: Distribution P(T) for AT^ = 4, a = 0.6, /3 = 0.2 and 
L — 400. The blue solid line is the whole distribution. The 
other curves are the contributions to P{T) corresponding to 
the passage of at most Nmax particles through the bottleneck 
during the time interval T. 

to it. We now consider the probability distribution P{T), 
shown in Fig. [6] 

We focus on the jammed state, though a similar anal- 
ysis could be done in the FF phase. If we assume as be- 
fore that a stationary state is established in the bottle- 
neck, then the probability for the bottleneck to be empty 
should be = at each time step and, if we neglect 
correlations between successive time steps, the distribu- 
tion P(T) = (1 - P0)^-iP0 follows. For = 1, the 
stationary state is established very rapidly. For larger A^5, 
non negligible corrections are present because there is a 
non vanishing relaxation time towards stationarity. The 
accuracy of our prediction for P{T) depends on the ratio 
between the typical times T and the relaxation time 0. 
The typical time T is expected to increase more rapidly 
with the bottleneck size A/'^ than the relaxation time 0, 
making our prediction more accurate for large A^^. In- 
deed, in Fig. [7|we observe that the relative error between 
our prediction and simulation results decreases again for 
large bottlenecks. For all A/'5, the prediction works better 
for small /3, and the difference between theory and numer- 
ics increases when the MC phase is approached. Indeed, 
when f3 increases, P^ increases, and the typical T value 
decreases. It can then become smaller than the relaxation 
time and a description of transients becomes necessary. 
Besides, the queues may not have enough time to form 
and thus do not overfeed the bottleneck. 

Another feature visible in Fig. [6] is that there is an im- 
portant contribution to P{T) from events during which 
only a single or very few particles pass through the bottle- 
neck. When these events occur, the system explores only 
very special configurations, which are not representative 
of the whole stationary distribution. Thus the system is 
still in a transient state when the next reversal occurs. 

As a conclusion the dynamics for the reversal of the 



p-4 



A bottleneck model for bidirectional transport controlled by fluctuations 




Fig. 7: Relative error between the numerical result and theo- 
retical prediction for the decay rate of the tail of P{T). Error 
is plotted against A^b, for various /3, L = 100 or 400, a — 0.6. 
Lines are guides to the eyes. 



explain the shape of the distribution around its maximum. 
A large part of the distribution is due to non-typical events 
where only a few particles go through the bottleneck. A 
complete understanding of P{T) should involve the study 
of non-stationary distributions and correlations between 
flux reversals. 

To conclude, the main feature that puts this new model 
apart from others in the literature is that fluctuations lo- 
calized in the bottleneck can have a macroscopic effect on 
the whole system. It provides a sensitive test for different 
theoretical approaches and can be easily tested numer- 
ically. While a phenomenological approach assuming a 
stationary state in the bottleneck gives surprisingly good 
predictions for the free flow and part of the jammed phase, 
some other observations trigger much more complex ques- 
tions involving non-stationary and correlated behaviors. 
While we concentrated on identical lanes without direc- 
tional bias, one can, of course, consider asymmetries either 
in particles species or capacities of the lanes which are of 
interest for different applications. 



flux actually involves two quite different mechanisms (one 
based on stationarity and the other involving transients), 
and exhibits some kind of intermittency, with long periods 
of one-directional flows alternating with rapid switches. In 
order to have a full description of the dynamics of the flux 
reversal, not only the transient nature of the incoming 
flows in the bottleneck should be taken into account, but 
also the correlations between successive switches. 
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Conclusions. — We have introduced a new model for 
bidirectional transport with a bottleneck. We have shown 
that theoretical considerations, based on the stationarity 
hypothesis for the flow inside the bottleneck, allow to pre- 
dict with a good accuracy the phase diagram, and the val- 
ues of the currents in the free flow phase, in the jammed 
phase, and in the MC phase for long bottlenecks. For 
short bottlenecks, transients dominate the behavior of the 
system, and as a consequence large values of the current 
can be observed, in spite of the cost of having to empty the 
bottleneck at each reversal of the flux. These transients 
could be studied by exploiting the exact results obtained 
for a single TASEP with a step initial condition [8)|20||2l]. 

The jammed phase turns out to have a striped struc- 
ture. In spite of the similarity with 19 , here the striped 
phase can be sustained in the bulk without any modifl- 
cation of the standard TASEP bulk rules. It should be 
noted that this striped phase results from a self-regulated 
dynamics for the flux reversal in the bottleneck. Though 
an assumption of constant reversal periods succeeds in pre- 
dicting the striped phase density, we flnd that actually the 
structure of the distribution of the bottleneck occupation 
times P{T) is more complex and emerges from two dif- 
ferent mechanism. The tail of the distribution can be ex- 
plained through our phenomenological approach assuming 
a 1-lane stationary state in the bottleneck. However, this 
prediction is valid only for small /3, and in any case cannot 
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